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D ' Abstract. The idea that the radio jets in AGN contain material in rel- 

ativistic motion is supported by many lines of observational evidence, 
\^ ' including morphology, brightness temperature estimated with interfer- 

ometers and with intrinsic variations, interstellar scintillations, X-rays, 
and superluminal motion. These are largely independent, and taken 
^ ' together make an irrefutable case for relativistic motion. 

00 ' 

in . 

\ 1 Introduction 

^ ' Scientists were surprised in the early 1950s, when it became clear that the "radio 
f^S ■ stars" in fact were extragalactic and hence very energetic, and the general idea of 
a calm orderly Universe had to be questioned. In later years the discovery of black 
£s . holes, quasars, pulsars, and other exotic objects changed the view completely, and 
brought us to the violent Universe as we understand it today. 

In a prescient article, Rees (1966) suggested bulk relativistic motion as an 
explanation for the rapid variations in flux density seen in quasars; and he also 
^3 ■ predicted the phenomenon of superluminal apparent transverse velocity, which was 
discovered five years later. At that point most people believed that cosmic jets 
contained relativistic jets of plasma, but what process could collimate and acceler- 
ate whole solar masses of material to relativistic speeds? The current paradigm is 
that a giant black hole lies at the center of a galaxy, and the process of accretion 
combined with a magnetic field generates the bulk relativistic flow. The details 
of this scenario are not well understood, but gravitational accretion is the only 
known means for generating the required energy per unit volume. 

In this paper we review the radio observations that provide evidence for rela- 
tivistic motion. The plan of the paper is to successively discuss frequency shift, 
morphology and Doppler boosting, brightness temperature measured with interfer- 
ometers and estimated with intrinsic variability, interstellar scintillations, inverse- 
Compton X-rays, and superluminal motion. Individually, and especially taken 
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together, these make a convincing case for relativistic motion in the radio jets. 
This conclusion is independent of any specific cosmology. 

Relativistic motion implies that the observed radiation patterns are narrow, 
that Doppler boosting increases the observed flux density, and that the observed 
time scale is shrunk below that in the frame of the AGN. These affect the ob- 
servable quantities; namely, frequency, luminosity, apparent transverse speed, and 
variability; and there should be statistical relations among them. We briefly touch 
on this in §4. 

A number of earlier papers have discussed the different types of observa- 
tional data bearing on relativistic motion. We mention only a few of them here: 



Ghisellini et al. (1993) compared X-ray and VLBI data to derive Doppler fac- 



tors and compared them with other beaming indicators; Vermeulen and Cohen 
(1994) assembled all the existing data on superluminal sources and analyzed sev- 
eral models involving the distributions of Lorentz factor and of the ratio of pat- 
tern to beam velocity; Lister and Marscher (1997) made Monte Carlo calcula- 
tions of distributions of the measureable parameters, and of their correlations; 
and Lahteenmaki & Valtaoja (1999) estimated Lorentz factors and also the in- 
trinsic temperature in the synchrotron sources. These past surveys mainly used 
inhomogeneous data culled from the literature. However, two large surveys on the 
VLBA, at 5 MHz |(Vermeulen et al. 2003)] and 15 GHz (Kellermann et al. 2004, 
hereafter K04) are now bearing fruit, and analyses with these superior data sets 
should soon give much better results. 



2 Formulae 

Here we list some formulae, for later convenience. Assume a relativistic luminous 
plasma blob with velocity (3 — v/c and Lorentz factor 7 = (1 — /3 2 )~2. If it is 
moving at angle 9 to the line-of-sight (LOS) then its Doppler factor is 



1 



7(1 — (3 cost 



(1) 



The observed frequency is v — 8v* , where the superscript * refers to the coordinate 
frame of the moving plasma. Non-starred quantities are in the frame comoving 
with the AGN. The brightness temperature has the same transformation law as 
the frequency: 

T b = ST b *. (2) 
The observed flux from a point source is 

s v = s 3 s;,, (3) 



while for a jet the exponent is 2, and in general the exponent can be between 2 
and 3[(Lind fc Blandford 1985] 
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If a source component has proper motion fi then it has an apparent transverse 
speed /3 app = /iD(l + z)/c, where D is the metric distance to the source. In terms 
of the source geometry, 



f3sin9 
1 — (3cos 

Graphs of this relation are in Figure |21 

3 The Observations 

3.1 Frequency Shift 

Several galactic sources show radio jets that are analogous to those seen in the 
high-luminosity extragalactic sources. SS 433 contains a jet with radio-bright 
components which move out from the core, and the jet changes orientation in a 
way that is consistent with precession. This source also exhibits frequency-shifted 
Balmer lines that, when tracked over time, show conclusively that the speed of the 
line-emitting clouds is f3 = 0.26 (Margon et al. 1980)| This value is not relativistic, 



but it is well above values otherwise associated with stars, including supernovae 
shock waves. If SS 433, a compact stellar object, can generate a jet with (3 = 0.26, 
then we might expect that an extragalactic object a billion times more luminous 
can produce a jet that is relativistic. 



3.2 Morphology and Doppler Boosting 

Early imaging showed that many radio galaxies contain lobes more-or-less sym- 
metrically located across the galaxy, and some objects also contain a bright point 
source at the center. A good deal of speculation went into the question of how 
these lobes were formed and energized. They were hundreds or even thousands 
of kiloparsecs from the center, but their synchrotron lifetimes were only a mil- 
lion years. If they were ejected from the galaxy it must have been at relativistic 
speed, otherwise there would be a lifetime problem. Alternately, there might be 
some kind of continuous beams that kept the lobes energized (eg Ryle and Longair 
1967). 

VLA images then showed that the beams did exist, and that there could be 
continuous feeding of the lobes. However, in most cases the beams were strongly 
asymmetric; one side typically was much stronger than the other. Another clue 
came from VLBI high- resolution studies. Many compact sources were linear but 
strongly non-symmetric; a bright "core" was at one end of a jet that had a spec- 
tral gradient such that the core had the flattest spectrum. See K04 for many 
examples. An important discovery was that the inner jet nearly always pointed in 
the direction of the stronger of the two outer jets. How could there be one-sided 
inner structure with two-sided outer structure? Most astronomers accepted the 
idea that the jets contained relativistic beams, and Doppler boosting enhanced 
the forward beam relative to the backward beam. 
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In 1987 the "Laing-Garington effect" was discovered (Garrington et al. 1988) 
this was essentially regarded as a proof of the relativity hypothesis. The two outer 
radio lobes of a radio source typically have different amounts of Faraday rotation, 
and the lobe with the weaker Faraday effect is nearly always on the jet side. It 
must be in front, because on average its radiation goes through less plasma than 
for the one in back, the non-jet lobe. The visible jet then is the one in front, and 
the only convincing reason that has been offered for this is Doppler boosting. 

The front-to-back ratio of the jets, r, is typically greater than 10. This does 
not call for a very high value of 7, although 3C 273 is exceptionally one-sided, 



with r > 5500 (Davis et al. 1985) The compact inner jets often have r = several 
thousand, which requires S > 7; i.e. 7 > 3.5 and 9 < 8.2°. However, these are the 
extreme values of 7 and 9, and more likely values are 7 ~ 7 and 9 ~ 6°. These 
calculations assume that the jet and counterjet have the same intrinsic luminosity. 

Note that the measured high values of r imply small values of 9. The strongly 
one-sided objects must be aimed close to the LOS. What of the mis-directed ones, 
those more numerous objects not close to the LOS? They are not Doppler-boosted; 
and in fact most will be deboosted and hence weak. The mis-directed objects 
are generally thought to be radio galaxies, and indeed these typically have weak 
cores (Urry & Padovani 1995) The strong, compact, one-sided, flat-spectrum 
sources often also show superluminal motion, which also requires an orientation 
close to the LOS. (see §3.6). Given the number of sources that are known to 
be highly aligned with the LOS, the parent population is likely to be very large 



(Lister & Marscher 1997) This topic, however, is beyond the scope of this review. 



3.3 Interferometric Brightness Temperature 

A relativistic plasma will produce synchrotron radiation at radio wavelengths 
and also X-rays, from the inverse-Compton process. For reasonable values of 
the magnetic field the maximum intrinsic temperature that the plasma can have, 
if it is in dynamic equilibrium, is T Q ~ 10 12 K, because at a higher tempera- 
ture second order scattering becomes important and rapidly quenchs the plasma 



(Kellermann & Pauliny-Toth 1969) Somewhat similar arguments have been made 
by Singal (1986) and by |Sincell fc Krolik (1994)| who derive 3 x 10 11 K and 2 x 10 11 



K, respectively, as upper limits. However, Readhead (1994) has argued that be- 
cause at 10 12 K the magnetic and kinetic energies in the plasma are very far out 
of balance, a lower equilibrium temperature, T oq ~ 5 x 10 10 K, is more likely to 
be found. 

Measurements made with VLBI show that the temperature can be higher than 
these values. Figure QJi shows the distribution of Tb measured with the VLBA 



at 15 GHz (Zensus et al. 2002) and Figure ^> shows the distribution measured 



with the space VLBI program VSOP at 5 MHz |(Hirabayashi et al. 20 00) For 
both the VLBA and VSOP results, some of the measured values are lower limits. 
The two distributions are similar and show that some sources have Tb an order of 
magnitude and more above the SSC and other limits and estimates for Tb described 
above. Theoretical ways of exceeding the limits have been proposed, including a 
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Log Brightness Temperature (Source Frame) 



Fig. 1. Distributions of interferometric brightness temperature (a) measured with the 
VLB A at 15 GHz, N=19|(Zensus et al. 2002)| and (b) measured with VSOP at 5 GHz, 
N=59|(Hirabayashi et al. 2000) | Note the different vertical scales in (a) and (b). 



non-equilibrium plasma (Slysh 1992), and coherent radiation processes (Benford 
and Lesch 1998). However, we prefer (with most others) to invoke Eq (2) and 
explain the discrepancy with relativistic motion. We expect that most sources 
will have T Q ~ 10 11 K in accord with the estimates of Readhcad and others, 
although we also think there might be a distribution of T Q with a small number of 
sources approaching 10 12 . In any event, it is clear that many sources must contain 
relativistic beams with 5 = 5 — 10 and more. 

Note that this discussion and the one based on morphology are independent 
and complementary. They both require Doppler factors up to ~ 10. 



3.4 Variability Brightness Temperature 
3.4.1 Intrinsic Variablity 

Many AGN show rapid variability, with time scale r as short as a few weeks. With 
the assumption that r is limited by the light-crossing time, and with the aid of 
a model, r can be turned into an angular diameter. The "variability brightness 
temperature" T var can then be found from the flux in the outburst. However, 
Eq (2) is not valid for converting these T var to intrinsic values, T Q , because the 
transformation of r introduces another factor of 5, and r enters twice in calculating 
the solid angle. Hence, in this case (2) is replaced with 



T var = <S 3 T , (5) 

where T Q is the intrinsic temperature in the plasma. With the assumption of a 
particular value for T Q , Eq (5) can be used to calculate <5 va r; the cube root makes 
the estimate of 5 val somewhat insensitive to the model used in deriving the angular 
size. Lahteenmaki & Valtaoja (1999) have used their flux monitoring data at 22 



and 37 GHz to calculate T var and 5 val , for T Q = 5 x 10 10 K. Their histogram of S v 
is in Figure [3 Their values have been converted to the cosmology we use in this 
paper, H Q = 70 km sec -1 Mpc -1 , fl m = 0.3, f^A = 0.7. The distribution of i5 var 
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Variability Doppler Factor 



Fig. 2. Distribution of Variability Doppler Factor from Lahteenmaki and Valtaoja (f999), 
with their cosmology corrected to the values used in this paper. N=81 

decreases steadily with 5 var up to <5 va r ~ 30. The model used in these calculations 
is an optically thick sphere of radius cr. 

If the interferometric data in Figure^are converted to 5 with the same intrinsic 
temperature as used for the variability, T Q = 5 x 10 10 K, the 8 are larger than 
the ones in FigureEl and a higher temperature, around T = 10 11 K, is needed to 
bring them into rough agreement. A detailed comparison cannot be made because 
the selection criteria for the two data sets are very different; and also because the 
interferometric data contain a number of lower limits on Tb- 



3.4.2 Interstellar Scintillations 

Some compact sources show rapid variability due to the passage of the Earth 
through the irregular diffraction pattern produced by the interstellar medium. 
This interstellar scintillation can be as fast as a few hours, giving intraday vari- 
ability, or IDV. Scattering theory gives an estimate of the distance to the screen, 
and the angular size of the background source. Thus the brightness temperature 
can be estimated; the method is analgous to VLBI, but with a longer baseline. 
Analysis of the IDV cases shows that the screen distance is only 10s of pc, and 



the angular scale is 10s of micro-arcsec (Dennett-Thorpe & de Bruyn 2000) The 
best-determined case appears to be PKS 0405-385 |(Rickett et al. 2002)[ "^hich has 
T b ~ 2 x 10 13 , and 6 ~ 75 (for T a = 3 x 10 u , the value used by Rickett et al. The 
value 6 ~ 75, while high, is not unprecedented, because the generally-accepted 
theory for 7-ray bursts requires jets with 8 ~ 100 or more. There are many AGN 
jets with 5 up to about 30, (Fig |2J) and a few with 5 near 100, but none have 



been found in-between. Lahteenmaki & Valtaoja (1999) state that, if there were 



sources with variability time scale r ~ 1 week at 22 and 37 GHz, they would have 
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Fig. 3. Apparent transverse velocity /3 apP vs 8, the angle to the line-of-sight. Note 
that every /3 app > 1 has a maximum possible angle, sin^max = 2/3 app /(l + 0i pp ) and a 
minimum Lorentz factor 7 m i n = (/3a PP + l) 1 ^ 2 



been found. However, in view of the selection biases in all these studies, we think 
it premature to generalize these statistics. 



3.5 Inverse-Compton X-rays 

A synchrotron source in dynamic equilibrium must radiate X-rays via the inverse- 
Compton process, but in some cases the measured X-rays are much lower than 
the value expected from the radio flux density. The accepted explanation for this 
is relativistic motion, and S can be estimated if the appropriate radio and X-ray 



measurements are made (Marscher 1983) The calculation is strongly dependent 



on the radio diameter, which has to be measured with VLBI. Thus this method 
is closely allied to the VLBI temperature method (§3.3). The X-rays should be 
measured simultaneously with the radio, and the radio observations should be 
made at the peak of the synchrotron spectrum. These conditions are seldom 
met. Nevertheless, calculations made from literature surveys (eg Ghisellini et 
al. 1993) show that many flat-spectrum sources appear to require 5^5 — 10. 
Bloom et al. (1999)[ however, show that the errors typically are fairly large, and 



only a few of the sources in their sample actually require S > 1. The best case 
appears to be 3C 345; Unwin et al. (1994, 1997) show that in this source 5 increases 
from ~ 5 to > 10 as the radiating component moves out from the core. (These 
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12 3 

Redshift 



Fig. 4. Proper motion vs redshift (fiz diagram) for the fastest "Excellent" or "Good" 
component for 110 sources in the 15-GHz VLBA survey. The line 7 = 25 is the locus of 
points where the minimum value of 7 is 25; a source on the line can have a higher, but 
not a lower, value of 7. 



values become somewhat larger if the the cosmology used in this review is used.) 



3.6 Superluminal Motion 

A local observer sees the radiation from a relativistic jet contracted in time, pro- 
vided S > 1, and the apparent transverse speed, /3 app , is greater than the speed of 
light, in some range of 6 (Eq. 4). Figure shows /3 app vs 9 for several values of 7; 
for 7 2 >> 1, /3 app is a maximum at 9 w 1/7, and /3 app , ma x ~ 7- Conversely, for a 
given /3 app , we have 7 min ~ /3 app . 

Early VLBI observations showed a number of sources that displayed this super- 



luminal effect (Vermeulen & Cohen 1994) More recently, systematic large-scale 



surveys have produced more than 100 well-defined superluminal sources. In this 
section we show some results from a survey on the VLBA, and then discuss alter- 
natives to the relativistic beam for explaining the observations. 

Figure 0] shows the nz diagram for the fastest components that are rated E 
(excellent) or G (good) for 110 sources in the 15-GHz VLBA survey (K04). The 
solid line, marked 7 = 25, shows the points (/J.,z) for which 7 > 25. Hence, the 
points close to the line have 7 rj 25 or greater, and sources below the line can have 
correspondingly lower values of 7. This plot shows directly that the measured 
patterns must be moving relativistically. Figure shows the histogram of /3 app for 
the sources in Figure 0] This also is a histogram for 7 m ;„. 

We note that the fastest component found in even a rather small sample 
of beamed sources with constant 7 will have /3 app which closely approaches 7 



(Vermeulen & Cohen 1994) Hence, we take 7 max = 25 in our later discussions. 



M.H. Cohen et al: Relativistic Phenomena in Active Galactic Nuclei 



60 




Fig. 5. Histogram of the apparent transverse velocity, /3 apP , for the 110 sources in Fig- 
ure 01 The median value of /3 app is 4.2. This also is a histogram of the minimum possible 
value of 7. 
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Fig. 6. Histogram of # max , the maximum possible angle to the LOS, for the 1 
Figure [SJthat have S > 1.5 Jy at 15 GHz. The peak near # max = 90° consists of sources 
with low or non-significant values of /3 app which do not constrain the angle. 



In reality the value could be somewhat larger, but not smaller. 

Any /3 app > 1 has a corresponding maximum possible 0, sin6* max = 2/3 app / (1 + 
fllpp) ( see Figure EJ. The histogram of 6> max for the 68 sources with S > 1.5 
Jy is in Figure HO We picked this flux density cutoff because it is the limit of 
the original selection, and the weaker sources probably introduce a serious bias. 
These 68 sources, roughly speaking, approximate a flux-limited sample. In plotting 
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this histogram, we set every source with /3 app < 2a or /3 app < 1 to 90°; this is 
responsible for the peak at 90°. 

In fact, 6* max is highly unlikely to be the actual value of 9 for a beamed source. 
It can be shown that the most probable viewing angle, 9 P , for a source picked on 
the basis of its flux density from a randomly-oriented beamed population, is about 



1/(27) (Lister & Marscher 1997) Since the minimum 7 is approximately /3 app , a 
rough upper limit to 9 P is 6* max /4. We conclude that the histogram in Figure is 
unlikely to be representative of the real values of 9; a more probable distribution 
would be more peaked at small angles. 

Lister and Marscher (1997) show some plots of the distribution of /3 app and 
9 from Monte-Carlo simulations of observing a complete well-defined sample of 
relativistically boosted jets. We have begun to observe such a sample, in a con- 



tinuation of the 15 GHz VLBA survey called the MOJAVE project (Lister 2003) 



3.6.1 Pattern Speed vs Beam Speed 

It is usually assumed that the speed measured for the moving pattern is also the 



speed of the beam, but this often may be wrong (Lind & Blandford 1985) Shocks 



will exist in the jet, and in some cases the moving blobs will be compressions 



due to the shocks (Hughes et al. 1985) In these cases the pattern speed can 



be different from the underlying flow velocity. As an extreme case, consider a 
standing shock wave set up at a bend in the jet; the measured speed (zero) gives 
no information on the actual flow velocity. Another possibility is that a beam 
contains two components, a narrow fast spine and a slower sheath. In this case 
the observed pattern speed might be a combination of the two speeds, although 
we might see only the spine, because its Doppler boosting will be higher. Other 
effects, including bending of the beam and evolution in the moving component, 
can also give erroneous speeds. Nevertheless, in most cases the measured pattern 
velocity should be similar to the actual beam speed, because the speeds are roughly 
compatible with other relativity indicators. (See §4.4). 

3.6.2 Alternatives to the Relativistic Beam Model 

Relativistic beams are generally accepted as the explanation for superluminal mo- 
tion, but alternatives have been suggested and are still occasionally referred to. 
These can be grouped into two main categories. 
1. Random 

Among the 110 sources of the 15- GHz survey, nearly all show expansion or zero 
velocity. The two sources with a statistically significant (2a) negative velocity can 
be reasonably explained with outward helical motion or with a bent jet, as an 
alternative to contraction. Further, in most cases the motions are along a line, 
and the expansions can persist for up to a decade in distance from the core. In 



the two-sided source NGC 1052 (Vermeulen 2003) the speeds are similar on the 
two sides. These motions are not random. 
2. Non-beaming models 
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The light-echo model (Lyndcn-Bell 1977) works for supernova remnants but 
not for jets in AGN, because the predicted distribution of velocities does not 
match the data. The same problem arises for echoes in a dipole magnetic field 



(Bahcall & Milgrom 1980) Ekers & Liang (1990) produced an ingenious model of 



a relativistic excitation wave in a jet, where the radiating material is stationary; 
i.e. it is not swept up by the wave. This suffers from the same problem with the 
distribution of velocities. 

These models do not have relativistic motion of the radiating material, so 
there is no Doppler favoritism. The selected sources must be oriented randomly 
relative to the LOS. On the other hand, the distribution in Figure® based only on 
the observed values of /3 app and without explicitly or implicitly assuming Doppler 
boosting, shows that the orientations are not random, which would give a sin# 
distribution. 

This result was quantified by Cohen (1990) with an analysis of the 12 strongest 
sources that were then known. Each source has a solid angle associated with max 
that can be compared with 2ir, the total available, in a manner analagous to a 
V/V m calculation. The result was that the probability that the objects were drawn 
from a randomly oriented parent population was small unless the Hubble constant 
was well above any reasonable value. This calculation can now be repeated with a 
much better and larger sample, and with a current value for H Q . The result for the 
fastest component that is rated E or G, for the 53 quasars that have S(15)> 1.5 
Jy (the original limit of the 15-GHz survey), is that the probability that they 
were drawn from a randomly oriented parent population is p < 6 x 10 -6 . This is a 
strongly biased estimate because (a) a distribution based on a more probable angle, 
rather than the maximum angle, would have a smaller probability, and (b) every 
source with /3 app < 2a or /3 app < 1 was set to 90 degrees. Thus the estimate for p 
is conservative, and we conclude that the quasars are not oriented at random, but 
rather are preferentially close to the LOS. The models without Doppler favoritism 
are not viable. 

An exception to this conclusion could arise if there were obscuration near the 
equatorial plane, so that small values of 9 would be favored. However, we do not 
see how this could come about. The dust clouds invoked for optical obscuration 
are transparent at radio wavelengths, and free-free absorption at 15 GHz would 
require unreasonable densities and masses of plasma. 



4 Self-Consistency 

The brightness temperature of a source can be determined with intcrfcrometry 
(Tb) and also with flux density variations (T var ). Comparisons of the two measures 
can help determine if the various assumptions underlying these determinations are 
valid. If Tb and T var are measured then the Doppler factor 5 and the intrinsic 
temperature (T Q ) can be calculated with Equations 2 and 5. Lahteenmaki et 
al. (1999a) have done this for several sets of sources, and found that most of the 
T Q have a range around 10 11 K. This is rather lower than the SSC limit but close 
to the other values discussed in §3.3. This suggests that, although the detailed 
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variability model commonly used is not too realistic, the results none-the-less 
should have some reliability, and the different measurements pertain to the same 
relativistic plasma. As discussed in §3.4.1, the sets of 5 found by the two methods 
are in rough agreement with T Q = 10 11 . 

Now we consider the relation between the variability Dopplcr factor S var and 
the apparent transverse speed /3 app , and assume that the same plasma motion 
affects both phenomena, through the shrinkage of the time scale by the forward 
relativistic motion. They should be closely connected, but note that high <5 var can 
go with low /3 app , when 9 < 1/7. If <5 var and /3 app are known, we can calculate 7 
and 9 from Eqs 1 and 4, but there is ambiguity because an assumed value of T Q is 
required for 5 var . We test the paradigm by calculating <5 var with several values of 
T , and comparing the resulting distribution of points (/3 app , S vaI ) with a Monte- 
Carlo simulation. This should show whether or not S var and /3 app could be related 
to each other through Eqs 1 and 4, and it might suggest an optimum value for T Q . 

The simulation is of the type described by Lister and Marscher (1997). It uses 
a set of beamed sources that are uniformly distributed in space and isotropically 
distributed in angle, and have power law distributions of luminosity and Lorentz 
factor. For the simulation shown here we use power law indices of -2.0 and -1.5, 
respectively. The calculation draws sources at random and counts those that are 
above a flux limit. 

Figure [3 shows /3 app plotted against S var for the 49 sources that have both an 
"excellent" or "good" valu e of /3 app (K04) and a value of <5 v ar from the 22 and 37 



GHz variability studies of Lahteenmaki & Valtaoja (1999) We have eliminated 



five sources with speed consistent with zero and located at a sharp bend in the jet, 
because we assume that they are connected to stationary shock waves and are not 
indicative of the motion. The Monte-Carlo simulation is in the top left panel; the 
other panels represent the observations with <5 var calculated with different values 
of T Q . These plots are discussed further in Cohen et al. (2003)] The solid lines 
in Figure [7| represent 7 = 25. A sample of sources with a distribution of 7 with 
7m ax = 25 will be entirely below the line. We used 7 = 25 following the discussion 
of Figure 0] and it appears to provide a reasonable envelope. 

The probability of finding a beamed source selected on the basis of flux density 
is maximized well inside the I/7 cone (Vermeulen & Cohen 1994) and 72% of the 
sources in the simulation are indeed inside the cone. In Figures Eb,c,d the fractions 
are 86%, 71%, and 63%, respectively; and on this basis we would select (c) as 
providing the best fit to the simulation. But the detailed distribution also suggests 
that (c) is preferable to (b) or (d), In (b) there is a cluster of points outside the 
7 = 25 circle including one at <5 va r = 70; however it is hard to judge the importance 
of this because the errors on <5 var are large and not well-understood. The effect 
becomes exaggerated if T D is reduced below 1 x 10 10 . In (d) the points appear to 
be clustered rather far to the left, and this is exaggerated if the temperature is 
raised. The result is that the speeds measured with the VLBA are consistent with 
the variability Doppler factors provided that the intrinsic temperature is roughly 
2 x 10 10 K. 

The number of points is small and the errors are large, but in Figure Et, 
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Fig. 7. Simulated and measured values of (/3 apP , <5 va r). (a) Simulated flux-limited sample 
(N=133) of beamed sources having 7 max = 25. The solid line is the locus of points for 
7 = 25; the dashed line /3 apP = 5 var marks the "I/7 cone" around the LOS (see text), 
(b) The points (/3 app , <5 var ) for the fastest component of the sources having both a /3 app 
value from the 15 GHz VLBA survey (K04), and a <5 var value from variability studies at 
22 and 37 GHz |Lahteenmaki fc Valtaoja (1999)| (N=49). The <5 var are calculated with 
intrinsic temperature T„ = 1 x 10 lu . (c) as in (b) with T = 2 x 10 10 . (d) as in (b) with 
T = 4 x 10 10 . 



for T Q = 2 x 10 10 K, the distribution is perhaps matched to the simulation. 
Cohen et al. (2003)| discuss the fit. T Q = 2 x 10 10 K can be taken as a rough 



value of the intrinsic temperature. In fact there surely is a distribution of T Q , but 
there are not enough data to study that. 

Figure |H1 shows the distribution of 7 calculated from the (/3 a pp , <^var ) points, 
with T„ = 2x 10 10 K. Four sources have formally negative values of /3 app , and in 
Figure |H1 they are plotted as upper limits at the 2cr position. The values above 
7 = 25 can be seen above the line in Figure [7] Two of them are within 2a of 
7 = 25, and the third, like all the points, has an unknown error connected with 
the variability Doppler factor. 

The value T Q ~ 2 x 10 10 K from the (/3 app , 5 var ) distribution is rather different 
from the value T Q ~ 10 11 K from the comparison of interferometer and variability 
Doppler factors. However, these agree within an order of magnitude and perhaps 
this is a success, given the many assumptions going into the analyses, and the 
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Fig. 8. Histogram of Lorentz factor calculated from <5 V ar and /3 apP with T = 2 x 10 K. 
Arrows indicate upper limits. 



large observational errors. 



5 Conclusions 

We have described five independent sets of radio phenomena which show that there 
is relativistic motion in the compact jets of AGN. These are the observations of 
morphology, interferometry, intrinsic variability, interstellar scintillations, and su- 
perluminal motion. X-ray measurements provide further confirmation, although 
these results are dependent on the VLBI measurements from which the brightness 
temperatures arc determined. Each of these phenomena has a theoretical diffi- 
culty which is readily alleviated by assuming that there is bulk relativistic motion 
towards the observer. The morphology, interferometry, interstellar scintillations, 
and the X-ray method are independent of redshift and hence of cosmology. Inter- 
pretations based on intrinsic variability and on internal proper motion, however, 
require a knowledge of the distance to the source. The variability methods, and 
the interferometry method, require an assumption of an intrinsic temperature in 
the synchrotron cloud. The distribution of S peaks at low values and decreases to 
5 ~ 30 The Lorentz factor 7 is a more fundamental quantity than 5, but is harder 
to estimate. Values estimated for 7 range up to 34. 

A comparison of /? app and <5 va .r shows that the superluminal speeds and the 
variability Doppler factors are roughly consistent with the standard relativistic 
motion paradigm, provided the intrinsic temperature in the synchrotron clouds is 
~2x 10 10 K. This value is below the various theoretical upper limits, and is close 
to Rcadhead's estimate (1994), based on equipartition between the energy in the 
particles and that in the magnetic field. 
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